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Abstract 

Pre-decoupling magnetic fields affect the scalar modes of the geometry and produce 
observable effects which can be constrained also through the use of current (as opposed to 
forthcoming) data stemming from the Cosmic Microwave Background observations. The 
dependence of the temperature and polarization angular power spectra upon the parameters 
of an ambient magnetic field is encoded in the scaling properties of a set of basic integrals 
whose derivation is simplified in the limit of small angular scales. The magnetically-induced 
distortions patterns of the relevant observables can be computed analytically by employing 
scaling considerations which are corroborated by numerical results. The parameter space of 
the magnetized Cosmic Microwave background anisotropies is also discussed in the light of 
the obtained analytical results. 
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1 Formulation of the problem 



There are two complementary approaches to the analysis of the Cosmic Microwave Back- 
ground (CMB in what follows) observables. The first one is direct and it consists in comput- 
ing the angular power spectra by faithfully including all the relevant physical effects. The 
second approach is indirect, i.e. it amounts to deriving the dependence of the (measured) 
temperature and polarization anisotropies upon the parameters of the the underlying model 
which needs to be falsified. The recent WMAP Syr data [H [21 [3] (see also [H [5]) have 
been confronted with a number of theoretical scenarios that are logically organized around 
the ACDM paradigm where A stands for the dark-energy component and CDM stands for 
Cold Dark Matter. Similar statements can be made for other recent CMB data such as the 
ACBAR observations [6l|7| and the QUAD measurements [H [9l lOl [11] . 

A useful bridge between the direct and the indirect approach is represented by a num- 
ber of scaling relations which serve as a diagnostic for the dependence of the (observed) 
angular power spectra upon the parameters of a pivotal model. The temperature and po- 
larization autocorrelations (i.e., respectively, TT and EE angular power spectra) and their 
mutual cross-correlations (i.e. the TE angular power spectra) can be written, with shorthand 
notation, as 

(TT) _ ^(^+ l) ^(TT) ^(EE) _ ^(^ + 1) ^(EE) ^(TE) _ (TE) . 

W - — ^ — ? - — ^ — > - — ^ — ■ 

In the ACDM scenario the angular power spectra of Eq. (11. ip are functions of, at least, six 
physical quantities 

= ^£ ^('^s, ^bO, ^cO;^A; Hq, t^e) , (1-2) 

where X and Y stand, respectively, for T and E and where the parameters denote, with 
standard notations, the spectral index of (adiabatic) curvature perturbations (i.e. rzg), the 
critical fractions of baryons, CDM and dark energy (i.e., respectively, f2bo, ^co and Qa), the 
Hubble constant Hq and the optical depth at reionization (i.e. ere)- 

In the ACDM paradigm as well as in it extensions, the known scaling relations are often 
not the result of a numerical inference but are derived by means of analytical methods. 
Suppose, for sake of concreteness, that all the parameters of Eq. (11. 2p are fixed to the best 
fit of the WMAP Syr data alone and just one (e.g. the spectral index) is allowed to scale. 
From semi-analytic considerations it follows that 

«r«(i)-\ ar^ijj^'\ cr^i^y ,.3, 

where the notation oc signifies that the corresponding quantity scales with the multipole 
in a given manneil^. When the scalar spectral index changes from the best-fit value (i.e. 



^In Eq. (II. 3|) £p denotes the pivot multipole at which the initial conditions are customarily set. This 
scale is largely conventional and it will be hereby chosen to coincide with £ — 29 which does correspond to 



the pivot wavenumber fcp — 0.002 Mpc ^. 
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rig = 0.963) to a different value G^^^^ and G^™^ will be modified according to Eq. (11.31) . On 
the vertical axis of the plots reported in Fig. [T]the ratios (ng = 0.963)/G^ {ris = 1) 
and Gf^Vs = 0.963)/Gf (rig = 1) are computed numerically (full line) and analytically 
(as they emerge from Eq. (II. 3p ). 

There are scaling relations involving, at once, different parameters. As it is known from 
elementary considerations, the height of the first peak in the acoustic oscillations of G^^'^^ 
scales with first power the sound speed of the baryon photon fluid, which depends, in turn, 
upon the critical fraction of baryons; in formulae: 

cgb(^.) = , ^ Rdz) = 7- = 30.36 (—), (1.4) 

y^3[l + i?b(^*)] \z,J 

where cUb = ^o^bo z^, is the redshift of the last scattering. The examples can be multiplied 
by considering all the parameters of the ACDM scenario either alone or in some appropriate 
combinations. Instead of considering the dependence of the temperature and polarization 
anisotropics upon the various parameters listed at the right hand side of Eq. (11.21) . it is often 
practical to consider a class of truly physical parameters emerging directly from the analysis 
of the various power spectra (see, e.g. [121 1131 [H [U] and references therein). Concrete 
examples along this line are: 

• the relative heights of the first three peaks in G'^^'^^] 

• the positions of the peaks in all the observed angular power spectra (i.e. G^^^\ G^™'' 
and G^'^^'*) their heights, their depths, their mutual distances; 

• the numerical value of the acoustic multipol^; 

and so on and so forth. The height of the first acoustic peak does not have a simple 
dependence upon the parameters of Eq. (II. 2p . Conversely, in terms of the quantities of 
Eq. (II. 2p . some power spectra exhibit rather contrived scaling properties which become 
instead manifest as a function of appropriate sets of derived variables which are accessible to 
direct observations. For instance, in the standard ACDM paradigm, the numerical values of 
the position of the anticorrelation peak in the TE power spectrum can be easily related to the 
position of the first Doppler peak in the TT power spectra; the height of the anticorrelation 
peak itself, however, does not have comparatively simple scaling with the parameters of Eq. 

The CMB observables can be indeed studied in terms of a set of so-called normal pa- 
rameters whose distinctive feature is that their mutual correlation is (or at least should be) 
very small. One of the purposes of the present paper is to look for similar types of scaling 
relations but in a qualitatively different case, i.e. when the model contains, on top of ACDM 



•^The acoustic multipole is defined as = pSz^,) Ir^iiz^) where DpSz^) and rs(z,t) are, respectively, the 
comoving angular diameter distance and is the sound horizon at last scattering. 
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Figure 1: The usefulness of scaling relations is illustrated in a simplified situation involving 
the variation of the scalar spectral index from the best-fit value of the WMAP Syr data alone 
(i.e. rig = 0.963) to an exact Harrison-Zeldovich spectrum (i.e. Ug = 1). 

parameters, also an ambient magnetic field. To be even more specific we wish to consider the 
situation where the ACDM paradigm includes also a magnetized background whose presence 
necessarily entails supplementary parameters. The minimal situation, in this respect, con- 
templates two new parameters, i.e. the magnetic spectral index ub and the magnetic field 
amplitude i?L- In this case Eq. f ll.2p becomes 

Cf^^^ = Gf^^\nB, Bi,, ris, J^bo, J^co,^a, -f^o,ero)- (1-5) 

To formulate in visual terms the main problem addressed in the present investigation, it is 
useful to look at Fig. [2] which, in some way, is the analog of Fig. [1] but in the case when 
large-scale magnetic fields are consistently included in the pre- decoupling physics. In the two 
plots at the left of Fig. [2] the TT and EE angular power spectra are presented in two cases, 
i.e. in the absence of an ambient magnetic field (dashed line in both plots) and in the case 
when a magnetic field modifies the initial conditions and the evolution of CMB anisotropics 
(full lines in both plots at the left). Just for illustrative purposes the magnetic spectral index 
and the comoving magnetic field amplitude have been chosen to be, respectively, n-Q = 1.5 
and -Bl = 20 nG. Always in Fig. [2] (but in the two plots at the right) the TT and the EE 
angular power spectra have been divided by the corresponding power spectra but computed 
in the absence of magnetic fields. In other words, in the plots at the right the two curves 
stem from the ratio of the angular power spectra illustrated in the right plots of the same 
Fig. [21 Already at a quahtative level, the right plot of Fig. [2] shows interesting features 
like, for instance, different periodicities but similar growth rates with the multipole number. 
Is it possible to understand these as well as other features in analytic terms? This is one 
of the questions we ought to address. Another question could be: as in Fig. [1] the scaling 
properties can be used to infer the variation of the power spectra with the scalar spectral 
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= 71 .9, a = 0.0441 , n „ =0.21 4, =0.742, e = 0.087 = 71 .9, a „ =0.0441 , Q. „ =0.21 4, =0.742, e = 0.087 

' bO ' cO A 're ' bO ' cO A 're 




Figure 2: The temperature and polarization autocorrelation (plots at the left). In the plots 
at the right we report the ratios Gf^^\nB 7^ 0, i?L 7^ with X = T, E. By definition 

^ixx) ^ g(xx)^^^ = 0, i?L = 0) i.e. Gf ^ denotes either the temperature or the polarization 
autocorrelations in the absence of ambient magnetic field. 

index, can we do the same when large-scale magnetic fields are present? To achieve such 
a purpose it will prove useful to employ different approximations schemes which served as 
basic ingredients for developing the numerical techniques that led to the results of Fig. [21 

Before going through the plan of the investigation it is appropriate to swiftly remind 
the main motivations related to the study of the magnetized Universe (see, e.g. |16l [17] 
for dedicated reviews on the subject) which is basically the research program pursued here. 
Magnetic fields in gravitationally bound systems are a rather mundane feature of our Uni- 
verse but one of the least understood especially when the typical correlation scale of the field 
is large as it happens, for instance, in the case of galaxies, clusters or even superclusters. 
Since the pioneering works of Alfven and Fermi fiSi [T9l [20] (see also [211 122] ) large-scale mag- 
netic fields have been the subject of numerous speculations whose detailed nature cannot be 
summarized here. The interested reader is refereed to review articles such as [TB] and to the 
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references of more recent publications [221 EH ESI EHl EZI ■ It would be difficult to account for 
the detailed arguments leading to the different ideas which can eventually originate large- 
scale magnetism. Furthermore, as it will be argued below, the overall spirit of the present 
approach is more pragmatic; instead of dwelling on the rigor of the speculations leading 
to large-scale magnetism, it seems more urgent, on a physical ground, to decide which are 
the measurable effects of pre-decoupling magnetism. The latter quest is also experimentally 
better founded since the degree of rigor of a speculation is always rather complicated to 
assess and might strongly depend upon the theoretical prejudice of the authors. 

The characteristics of the approaches discussed in this paper do not exclude the possibility 
of magnetizing the early history of the Universe [TB] which is a rather intriguing subject 
of speculations dating back to the pioneering works of Zeldovich [2BJ and, independently, 
Harrison [291 SO] (see also [16] for further details). The early phases of the evolution of 
the Universe are often connected with an inflationary epoch where, however, it seems 
to be rather difficult to produce large-scale magnetic fields in four- dimensional Friedmann- 
Robertson- Walker Universes. This impasse is related to the peculiar form of the evolution 
equations of Maxwell fields in curved backgrounds which are, technically, invariant under 
the Weyl rescaling of the geometry. To amplify gauge fields one might want to extend his 
model to higher-dimensional frameworks [31j or couple the kinetic term of the gauge fields 
to a spectator field |23j. Magnetogenesis models based on the dynamics of an appropriate 
spectator field delicately improve on the structure of conventional infiationary models and 
can be directly constrained by CMB data [23]. In spite of the specific model it seems relevant 
that the amplified gauge fields are Abelian. The only non-screened vector modes that are 
present at finite conductivity are the ones associated with the hypercharge field. The non- 
Abelian fields develop actually a mass and they are screened as the Universe thermalizes. 
After the electroweak phase transition the photon field remains unscreened with amplitude 
cos 9^y where y is the hypercharge field and 6'„ is the Weinberg's angle. While it is certainly 
interesting to speculate on the origin of large-scale magnetism prior to matter-radiation 
equality (i.e. for redshifts larger than, approximately, 3200) it is also rather urgent, as it 
will be argued in section 2, to scrutinize the CMB observables with the purpose of falsifying 
the statement that large-scale magnetic fields were indeed present around matter-radiation 
equality and, later, at the epoch of photon decoupling (i.e. for a typical redshift of the order 
of 1100). 

Large-scale magnetic fields are a well defined object of experimental study since the 
pioneering contributions of Hiltner and Hall correctly interpreted in terms of a large-scale 

Large-scale magnetic fields produced inside the Hubble radius after inflation will have a correlation 
scale bounded (from above) by the Hubble radius at the moment when some charge separation is produced 
(be it, for instance, the electroweak time). Since the Hubble radius, during radiation, evolves much faster 
than the correlation scale of the produced field, the typical scale over which the magnetic field is coherent 
today is much shorter than the Mpc, obliterating, in this way, the possibility of successfully reproducing the 
galactic magnetic field [16]. 
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(interstellar) magnetic field by Davis and Greenstein [32] • For extended reviews on galactic 
magnetism see [33] and [31]. Large-scale magnetism is also a well established phenomenon 
in rich (i.e. Abell) clusters of galaxies [35] (see also [361 EZ]) cis well as, still with large 
uncertainties, in superclusters [38]. Magnetic fields in normal galaxies at high redshifts [39] 
could be already strong and this would be theoretically suggestive, in spite of the necessary 
caveats stemming from the large observational uncertainties (similar, in nature, to the ones 
experience while observing supercluster magnetism). Since we do observe large-scale for 
progressively larger redshifts it is natural to ask what happens at the photon decoupling 
especially because of the wealth of CMB data. In investigating such a class of phenomena 
the idea is to frame the least number of hypotheses on the subsequent evolution of large-scale 
magnetic fields so that we will take as starting point a faithful plasma description which is 
exactly the one employed in terrestrial laboratories [101 HH HJl HS] . 

Having spelled out the general perspective of the present paper, its layout is, in short, 
the following. In section 2 the physics of the strongly and weakly interacting species at 
the epoch of photon decoupling will be briefly summarized with emphasis on the role of 
large-scale magnetic fields. Such a discussion will provide, in a reasonably self-contained 
perspective, all the equations which will be employed in the subsequent analysis. Section 
3 treats the analytical methods employed in the line of sight solutions of the Boltzmann 
hierarchy, while, in section 4 the (magnetized) temperature and polarization anisotropics 
are computed in terms of a set of 8 basic integrals. Section 5 illustrates the comparison of 
analytical and numerical results while Section 6 deals with the analysis of the parameter 
space of magnetized CMB anisotropics. Section 7 contains the concluding considerations. 

2 Strongly and weakly interacting species 

Prior to decoupling the evolution of the plasma can be described in terms of charged and 
neutral species. Charged species interact directly with the ambient magnetic field whose 
effect, on the neutral species, is mediated by the relativistic fluctuations of the geometry. The 
separate role of the electrons and ions is often overlooked when the relativistic fluctuations 
of the geometry are consistently taken into account. The approach initially formulated in 
[1^ US] and developed in [IH] is based on a rather conservative perspective: instead to 
doing a lot of effort to compute exotic phenomena triggered by large-scale magnetic fields, 
it is logically more urgent to compute in detail how large-scale magnetic fields affect CMB 
observables (see [17] for a more detailed formulation of such a research program). If large- 
scale magnetic fields gravitate and interact, simultaneously, with electrons and protons the 
most relevant effects on the temperature and polarization observables will be the one coming 
from the scalar modes of the geometry [HI H^l 1^ . 

The scalar modes of the geometry admit two kinds of initial data which allow for the 
inclusion of large-scale magnetic fields, i.e. adiabatic initial conditions leading to the mag- 
netized adiabatic mode (see [HI HHl HH]) and entropic initial conditions leading to various 
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magnetized isocurvature modes (see [H]). In what follows the main focus will be on the 
magnetized adiabatic mode. However, the same considerations developed here can be easily 
extended to the case of entropic initial conditions. In the case of adiabatic initial conditions 
the fluctuations of the spatial curvature are due to the fluctuations of the energy density 
while the fluctuations of the specific entropy are strictly vanishing at large scales. In the 
case of entropic initial conditions the curvature inhomogeneities are due to the fluctuations 
of the sound speed which are related, in turn, to the fluctuations of the specific entropy (see 
last part of section [2]). Non-adiabatic initial conditions can be observationally constrained in 
different ways [501 EH [521 EH [54] and may lead, in the magnetized case, to interesting shape 
effects on the CMB observables. There are, of course, also different themes which involve 
the physics of large-scale magnetic fields in connection with CMB physics (see [17] for a ded- 
icated review). For instance, large-scale (tangled) magnetic fields might have also specific 
effects related to the vector and tensor modes of the geometry (which are minute at large 
scales). These effects have been analyzed, at various levels of concreteness, in [551 [SSI [57] 
(see also [SSI \5U\ and compare them to [III [i^ ligl [i^l [50]). 

Depending upon their interaction rates, the constituents of the plasma can be classified 
into two groups: the strongly interacting species (such as the electrons, the protons and 
the photons) and the weakly interacting constituents (such as the CDM particles and the 
neutrinos). The difference between the two aforementioned categories resides in electro- 
magnetic interaction which strongly affects the evolution of the electron-photon-ion system 
while it affects only indirectly the evolution of the weakly interacting species. The indirect 
effect of large-scale magnetic fields on the weakly interacting species comes from gravita- 
tional interactions: since large-scale magnetic fields gravitate, the relativistic fluctuations of 
the geometry are modified by their presence via the enforcement of the various constraints 
stemming from Einstein equations. There is a whole class of effects which are related to the 
high-frequncy branch of the spectrum of plasma excitations [60l [61] which can be treated 
within the same framework described here (i.e. the magnetized adiabatic mode). In what 
follows, however, the focus will be on the scaling properties of the TT, EE and TE correla- 
tions since the angular power spectra of the B-mode polarization have been the subject of a 
separate study (see [SO], second and third papers). It should be stressed that the values of 
the magnetic fields used in this paper are, sometimes, extreme, e.g. intensities of 10 nG are 
by far excluded both by direct limits stemming from the polarization observables [60l [61] 
and from the analysis of the peak structure of the TT correlations [62]. In [62] the WMAP 
5yr data have been analyzed by including the effects of large-scale magnetic fields. In this 
perspective, for instance, the parameters reported in Fig. [2] are excluded. More specifically, 
the valued (ne, -Bl) = (2, lOnG) are excluded, by the analysis of the TT and TE, to 95 
% confidence level. At the same time, it is useful to illustrate the results in terms of these 

'"^The magnetic field intensity and the magnetic spectral index are assigned as in |46l|48] (see also second 
and third paper of [60 ). In the present context, Bl is the comoving amplitude of the field regularized over 
a typical scale — Mpc~^. 
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extreme values since, in this way, the visual impact is more pronounced and the scaling of 
the results with the parameters of the ambient magnetic field more evident. 



2.1 Generalities 

The simplest description of the pre-decoupling plasma in the presence of large-scale magnetic 
fields can be derived from the general pair of equations: 

r;-\s';r = 87cGt;:, (2.1) 

V^F^" = Anf. (2.2) 

In Eq. (12.11) is the Ricci tensor, R is the Ricci scalar and is the total energy- 
momentum tensor of the system. In Eq. (12. 2p F'^'^ is the Maxwell field strength and j'^ is the 
total current of the system. Both the total energy momentum tensor and the total current 
must be covariantly conserved, i.e. 

V^j^ = 0, V.T^'' = 0. (2.3) 



In Eq. (12.31) denotes the covariant derivative. The total energy-momentum tensor is 
given by: 

= Ti;-^ + Ti;; + rit] + t^-^ + ri;^ + t^/^ + ri;^^^, (2.4) 

where the subscripts denote, respectively, the contributions of electrons, ions, neutrinos, 
photons and CDM particles. More quantitatively the energy-momentum tensors of the 
different species are: 

T"^ = P^^t)^U^ = PiM(l)Mfi), TJJ^ = PcM"c)wfc) (2-5) 

7-"/3 _ 4 (/3) _ a/3 rj^ap _ 4 (,3) p , . 



t; 



1 



(EM) - 4^ 



_p.,pP ^ l^^^pp^^p^j^ ^ pap ^ ^^^ap^ (2.7) 



where F^i = —a^ Si and F^ = —o^eijkB^ are the components of the electromagnetic field 
strengths expressed, respectively, in terms of the electric and magnetic fields. In Eq. (12.61) 
the energy-momentum tensor of the neutrinos should also contain a contribution from the 
anisotropic stress which is, however, fully inhomogeneous and affects the evolution of the 
curvature perturbations rather than the evolution of the background metric. The total 
current of the system is due to electrons and ions, i.e. 

9f^u Wfi) U"^) = 1, 9f^u wfe) = 1, (2-8) 
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where e denotes the electric charg^ ; fie and fii are the physical (as opposed to comoving) 
concentrations of the electrons and of the ions. 

The evolution equations of the background geometry follow directly from Eq. (12. ip by 
recalling that, in the ACDM paradigm, the geometry is conformally flat (i.e. g^^^, = o?ri^y 
where rj^j^y is the Minkowski metric): 

= ^aVt, (2.9) 
-n' = 47iGa\pt + pt), (2.10) 
p[ + 3nipt+Pt) = 0. (2.11) 

In Eqs. ([21]), (|2J0D and (12111) 

• the prime denotes a derivation with respect to the conformal time coordinate r; 

• Ti. = a'/a which also implies Ti. = aH where H = a/a (where the overdot denotes a 
derivation with respect to the cosmic time coordinate t; recall that dt = a(r) dr); 

• finally the total energy density and the total pressure are: 

Pt = Pe + Pi + Pj + Pu + Pc + Pa, (2.12) 

Pt = y + y-PA- (2-13) 

For purposes of presentation we started directly from the covariantly conserved evolution of 
the energy-momentum tensor. It can be shown that this description is fully equivalent to a 
truncated Vlasov-Landau description [iTj H8]. 

2.2 Strongly interacting species 

The photons, the electrons and the ions interact electromagnetically and their velocities are 
tied together by the presence of scattering terms. At the same time photons, electrons and 
ions affect the evolution of the background geoemetry (i.e. Eqs. (I2.9p - (l2.1ip and (I2.12p - 
( I2.13P ) and of its relativistic inhomogeneities. From Eq. (12. 2p the evolution of the Maxwell 
fields obeys 

V ■ E = A7ce{ni-ne), (2.14) 
V-5 = 0, (2.15) 
VxE=-B', (2.16) 

V X B = ATTe{niVi-ncVc) + E' , (2.17) 



^In this paper the units will be such that e^/ (he) = 1/137. Furthermore, as it is apparent from Eq. ()2.2|) . 
in front of Fa/iF"^, in the action, there is a factor 1/(16 tt) which is reflected in the An of Eq. (|2.2p . Within 
these conventions and imposing the natural system of units h = c = 1, 1 Gauss — 6.9241 x 10~^° GeV^. 
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where the comoving concentrations and the comoving electromagnetic fields are: 

rii = a^fii, Ue = a^Ue, E = a^£, B = a^B. (2.18) 

The evolution equations of the electrons, of the ions and of the photons must include the 
relevant scattering terms governing their mutual momentum exchanges: 

+ nv^ = —[E + t/e X 5] - V0 + ^^ar^eK - ^e) + aT,,{v, - tJe), (2.19) 

me a 6 Pe 

v! + nv, = -^[E + ViXB]-V^+ ^^aT^iiv^ - v^) + aT,,^{v, - v^), (2.20) 
nip a 3 Pi Pi 



1 



V. 



7 4 



-V5^ — V0 + ar^i('i7i — -ily) + aT^f,{ve — v^). (2.21) 



The relevant interaction rates between the different species appearing in Eqs. f l2.19p . (12.201) 
and (I2.2ip are given by 



e2\2 8 /e2^2 



= 'T-eCre7, T^i = '^i0"i7, (^e^ = -T^[ ) , '^i7=o^( ) ' (2.22) 

rei = ney^aei = rie, a.. = ^^lnAc, Ac = |^^, (2.23) 

where T is the temperature and Ac is the Coulomb log |10l HI]. In Eqs. (12.190 . (I2.20p and 
(I2.2ip denotes one of the two longitudinal fluctuations of the geometry [63j whose exphcit 
form is given by0 

4 ^00 = 2a2 0, 5,g,j = 2a^^5ij. (2.24) 
The density contrasts of the strongly interacting species evolve, respectively, as 

5' = -V ■ {/e + 3^' - -^E ■ We, (2.25) 
S'. = -V • tJi + 3^' + ^E ■ tJi, (2.26) 

S'^ = Ai;' - ■ V,. (2.27) 

Equations (l219ll . (Km and (K2T\f describe, together with Eqs. (12:251) -( 12:271) . a three-fluid 
system formed by photons, electrons and ions. At early times, i.e. well before photon de- 
coupling, Fei ^ ^ H and F^c ^ (since nip ^ m^). The three fluid system can 
therefore be described in terms of two effective fluids. The main equations of the system, 
in this regime, are therefore the appropriate generalization of the familiar magnetohydrody- 
namical reduction [lOl HH HJl HS] but in the case where the relativistic fluctuations of the 
geometry are consistently included in the original equations of the multicomponent plasma. 



''In Eq. (|2.24p and in what follows, Ss{...) denotes the scalar fluctuation of the corresponding quantity. 
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By summing and subtracting Eqs. fl2.19p and fl2.20p the following pair of equations can be 
easily obtained, i.e. 



J -X B - , 4 

— V0+-^< 

a^Ph 3 pb 

J = a{E + VkyiB), (2.29) 



vl, + nv,, = — V0+-^ar^e(u,-tTb), (2.28) 



where pb = {rric + m^)nQ. To derive Eqs. (12.281) and (12.291) it should be borne in mind that 
the plasma is globally neutral, i.e. that fii = n^. From the evolution equations of the density 
contrasts (i.e. Eqs. (12.251) and (I2.26P ) it follows that 



5; = -V-4 + 3^/^' + ^, (2.30) 



a^Pb 



where, by definition. 



5b = — — — 8e H — — ^i, 5b = — , 5pb = 5pc + 5pi- (2.31) 

mp + me m-p + rrie pb 

Equations ^M), (12:291) and fl230D together with Eqs. fl2:2T|) and ^Q?TJ\ describe the baryon- 
photon fluid whose velocities (see Eqs. (I2.2ip and (12.280 ) obey 

+ = - V0 + ^(t;^ - iTb), (2.32) 

0. Pb -n^b 

% = -^V5^-V(l) + e'{v,-v^), (2.33) 

where i?b is the baryon-to-photon ratio already introduced in Eq. (I1.4p and where the 
differential optical depth e' = aF-ye = ahoXede'y has been introduced. 

The evolution equations for the photon-baryon system are the basis for the magnetohy- 
drodynamical description of the problem and for the analysis of the initial conditions of the 
Einstein-Boltzmann hierarchy [44l l45l l46l |48] . The differential optical depth enters directly 
the visibility function which gives the probability that a photon is emitted between r and 
T + dr: 

/C(r) = e'e-<^'^'\ e(r, tq) = a{T')n,{T')a,jT' . (2.34) 

The visibility function which will be adopted for the analytic estimates can be approximated 
with a double Gaussian whose first peak arises around last scattering (i.e. for r ^ r*) 

/C(r) = J\f{a^) e , < r < r^, (2.35) 

where is an intermediate conformal time such that < where is the reionization 
time. In Eq. (12.350 . N{a^) is determined by requiring that the integral of /C(r) over r is 
normalized to 1. The WMAP data suggest a thickness (in redshift space) Az* ~ 195±2 which 
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would imply that a^,, in units of the (comoving) angular diameter distance to recombination, 
can be estimated as cr^/ro ~ 1.43 x 10~^. When tq ^ t^, and tq ^ a* the normalization 
appearing in Eq. (I2.35P can be estimated as N{a^) — > cr^T^ V^^/V. The second peak occurs 
for the reionization epoch. Also in this case the visibility function can be approximated with 
a Gaussian profile centered, this time, around t^c- The specific form of the profile can be 
obtained from Eq. fl2.35p by replacing (r^, a^) with {r^eci o"rec) and by taking into account 
that 2;i.e = 11 ± 1.4. The Gaussian (or double Gaussian) parametrization of the visibility has 
been used in several works (see, e.g. |M1 |65j and also [66| 167]). 

Prior to decoupling the system can be further simplified. The photon and the baryon 
velocities are quickly synchronized because of the hierarchy between the scattering rate and 
the Hubble rate. Thus, the evolution equations of the photon-baryon system effectively 
reduce to: 

5; = 4^'-^V-T;^b, (2.36) 

= 3^' - V ■ tJ^b + (2.37) 
a Pb 



i?b + l"""' P,(l + i?b) 4(1 + R^) + 4aV, 



In Eq. (12.381) . the shear viscosity term 77 = (4/15)p^Ath depends upon the photon mean free 
path Ath which is, in turn, inversely proportional to the differential optical depth. 

2.3 Weakly interacting constituents 

The effect of the ambient magnetic field on the weakly interacting species is mediated by 
the relativistic fluctuations of the geometry which are affected by the scalar modes of the 
electromagnetic background. The evolution of the CDM is given by 

5; = 3^' - V ■ tJc, K + = -V0. (2.39) 
The evolution equation of the neutrinos can be written instead as 

K = -\v-% + Ai^', (2.40) 

< = Va, - iv5. - V0, (2.41) 

< = ■ V,. (2.42) 

15 

In Eqs. f l2.4ip and (12.420 is related to the neutrino anisotropic stress as did^YV^ = (p^ + 
Py)V^ai,. The weakly interacting species are affected by the action of large-scale magnetic 
fields through the evolution of the fluctuations of the geometry which obey the (perturbed) 
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Einstein equations. The (00) and (Oi) perturbed components of Eq. (12.11) are, in the gauge 



Ex B] 



{Pt + Pt)vt 



0. 



(2.43) 
(2.44) 



The (ij) component of the perturbed Einstein equations can be broken, respectively, into a 
trace full and a trace less part: 

r + n{(j)' + 2^') + {n^ + 2n')(j) + - V) = ^irCa^iSsPt + SsPb + 4pe],(2.45) 



i,)- 



(2.46) 



where 6sPt and 6sPt are the total fluctuations of the energy density and of the pressure while 



SsPb 



n 



Svra^ 
1 



SsPe 



Svra^ 



SsPb 



SspB 



SsPe 



5sPe 
3 ' 



3 



n 



1 



47ra4 



3 



(2.47) 
(2.48) 



where, as in Eq. (12.181) . E and B denote the electromagnetic fields and, by definition, 
B = a^B and B^ = BiB^ while E'^ = EiE^. Furthermore, following the same notation 
employed for the neutrino anisotropic stress we shall denote 



didjll'i = {p^ + p^) VVe, didjlV^ = {p^ + p^) VVb 



(2.49) 



2.4 The magnetized adiabatic mode 

In what follows the attention will be focussed on the magnetized adiabatic mode for which 
all the possible entropic fiuctuations vanish. This requirement implies that 5^ = — 3(Ci — Cj) 
where the indices run over all the constituents of the plasma and where d = —if) + 5i/[3(wi + 
1)]. The vanishing of the entropy fiuctuations vanish at large-scale (i. e. for kr <^ 1) implies 
that (i = Cj, where, again the indices run over all the species of the plasma. The latter 
(gauge-invariant) condition reads off, in the gauge defined by Eq. (I2.24p as 



di = —d 

+ 1 



J' 



(2.50) 



for any pair of constituents of the plasma. 

To set the initial conditions of the Einstein-Boltzmann hierarchy the consistent solution 
of Eqs. (1236|) -( 12:38D . (1239D . (ICTD -f irail . ([2331) -(EH]) and (IXiBD - dMD should be found 
at the initial integration time Tinitiai when the wavenumbers satisfy A^Tinitiai ^1. As it is well 
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known the latter condition implies that the corresponding wavelengths are larger than the 
Hubble radius: 

in Eq. (12.511) . k is the comoving wavenumber and /cphys(T) = k/a{T) is the physical wavenum- 
ber; furthermore 7i = aH where H = a/a and Ti = a' /a (see the comments after Eqs. (12. 9p . 
(I2.10p and (12.111) ). Equation (I2.5ip stipulates that the physical wavenumbers, at a given 
time, are always smaller than the Hubble rate implying, by definition, that the correspond- 
ing (physical) wavelengths are larger than the Hubble radius H~^. Note that the first equality 
in Eq. (I2.5ip is exact in a pure radiation-dominated phase when Ti = r^^. In the realistic 
situation, however, the scale factor interpolates between the radiation-dominated and the 
matter-dominated epochs and Ti = {2/ri)\/a + 1/a where a = a/a^q = [{t/tiY + 2(r/ri)] 
(see also the discussion around Eqs. (13. lip and (13.121) ). 

If the condition (I2.5ip holds, then it is also true that k/a <^ k/Ti where a = a^a is the 
conductivity. Indeed recall that 



^ I ^ ^ H = J^^-^, (2.52) 



ttem \meaj V 45 Mp 

where T = aTphys, gp is the effective number of relativistic degrees of freedom and Mp 
1/VG ~ 1.22 X 10^^ GeV. Then we can write, in the case of a cold plasma, 

k ^phys / fT^e / ^phys \ / 47r'^(y'p \ / 171^ /Tphys 



a "'^'"TpCs V Tphys - "^"^l i/ j I 45 j UphysJ I Mp j- ^^-^^^ 

According to Eq. fl2753|) . for Tphys < MeV, k/a < 1 provided kp^^jH < 1. Equation fl2753D 
shows also that that the condition k/a < 1 holds also for wavelengths which are shorter than 
the Hubble (i.e. /cphys > H) since Tphys/Mp is a really minute number (of the order of 10~^^ 
for temperatures in the eV range). 

The approximation scheme defined by Eqs. (I2.5ip . (12.521) and (12.531) allows for a con- 
sistent solution of Eqs. fl236D - (12381) . (I239|) . fl2:40|) - (l2:42|) . I^M)-(^M) and ([225])-([236D; 
the explicit form of the magnetized adiabatic mode can then be written as: 

' 4i?^ + 15 4i?^ + 15 

V'(fc) = (l + '^Ru^m + ^{i[Mk) + Mk)] - Rui^Bik) + QEik)]}, 
6^{k, t) = -20(/c) - i?^[OB(fc) + n^ik)], 
5,{k) = -2(l){k) - R^[nB{k) + fiE(fc)], 

Uk) = - ^^7[^b(A:) + QEik)], 

W = ~m - 7^7[^b(A:) + nEik)], 
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Mk, r) = -^[a^{k) + aE(fc)] + - mi 



euik,T) 



2 



m-^'^ + ^fiaM+Mk)) 



e,{k,T) = '^m, (2.54) 
where TZ^,{k) is the curvature perturbation on comoving orthogonal hypersurfaces, 

T^Ak) = - ^ -m - (2.55) 

For notational convenience, in Eq. fl2.54p . = 1 — Ry denotes the photon fraction in the 
radiation plasma and Ry is given, by definition, as 

p^ + py 1 + 3 X (7/8) X (4/11)4/3 ' ^ ^ 

where 3 counts the degrees of freedom associated with the massless neutrino families, (7/8) 
arises because neutrinos follow the Fermi-Dirac statistics; the factor (4/1 1)^/^ stems from 
the relative reduction of the neutrino (kinetic) temperature (in comparison with the photon 
temperature) after weak interactions fall out of thermal equilibrium. 

In Eq. (12.541) the following dimensionless quantities have also been introduced: 

He = = (2.57) 

see also Eqs. (12.471) for a definition of 6psPE and 6sPb', the quantities cte and ctb have been 
already introduced in Eq. (12.491) . It is finally useful to recall a pair of useful vector identities 
which connect cte and ctb to ^e and Qb, i-e. 

VVe = + ^J V ■ [(V xE)xE]- AttE ■ VpJ, (2.58) 

VVb = ^ + ^V-[Jx5], (2.59) 

where Pq = e{ni — ric); J = V x B/^An) is the total current in the one-fluid limit; R^ (see 
Eq. (II. 4p ) is the baryon-to-photon ratio. 

Because of Eqs. (I2.52p - (l2.53p . the contribution of the electric field fiuctuations to the 
initial conditions turns out to be almost always negligible. At the same time it is interesting 
to consider, in some detail, the transient regime where some putative electric field dies off 
thanks to the large values of the conductivity. In the latter case Eqs. fl2.54p and (12.551) 
allow for the inclusion of the electric field spectra in the initial conditions of the Einstein- 
Boltzmann hierarchy. This possibility will not be considered here but will be separately 
discussed. 
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2.5 Line of sight solution of the Boltzmann hierarchy 

The temperature and polarization power spectra are, by definition, 



r^(TT) 

— 


1 


m 


2£+ 1 


^(EE) 

- 


1 


m 


2i+l 


^(TE) _ 


1 


m 


2i + l 



(2.60) 
(2.61) 
(2.62) 



In terms of tlie intensity and polarization fluctuations in real space (i.e. Ai(n,r) and 
AE(n, r)), the coefficients a^^. and af^^ are: 



(T) 



dhY;^{n)Ai{h,T) 



(2.63) 



where n denotes the direction of propagation of the radiation and A*"^ = y'(£ — 2)!/(£ + 2)!. 
The (real space) E-mode fluctuation is defined as [TQl [71] (see also [60t 172]): 



AE(n,r) = -i{Ki^)(n)[Ki')(n)A+(n,r)] +i^i-^)(n)[i^i-')(n)A_(n,r)]}, (2.64) 

where A±{h, r) = AQ{n, r) ± z Au(?^, t) and where K^{n) are a pair of differential operators 
which can either raise or lower the spin-weight of a given function: 



±s 



9,9 ± 



smv 



(sin-i?) 



(2.65) 



The known advantage of dealing directly with the E-mode polarization (rather than with 
the Stokes parameters) is that AE(n, r) is a scalar (i.e. a function of spin-weight [71]) for 
rotations around the direction of propagation of the radiation field. In this sense Ae (n,r) 
is fully analog to Ai(n, r), i.e. the brightness perturbation of the intensity of the radiation 
field which is, of course, a function of spin-weight 0. Defining the projection of the Fourier 
mode in the direction of the photon momentum as /i = cos i), the evolution of the intensity 
of the radiation field reads (in Fourier space): 



a; + (ikfx + e')Ai = Si{k, fx, r). 



where 



Sp{k,T) = Ai2 + Apo + Ap2 

Si{k, fi, t) = —ikfi(f) + ip' + ^' 
e' 

Svd^, r) = - — (3«Aii + Vb) 



(3/^2 - 1) ^ , 
Aio + Ait-b bp{k, T) 



4 



ik 



(2.66) 
(2.67) 

(2.68) 
(2.69) 

(2.70) 



17 



Equation fl2.67p is the Fourier space version of Eq. (I2.32p and vy, is the divergence-full part 
Ub- The Fourier transform of the intensity and of the E-mode polarization is defined, within 
the present conventions, as 

Ai(ri, r) = J £k Ai(fc, /x, r), AE(n, r) = J d^k AE(fc, /x, r). (2.71) 

From Eqs. ( I2.64p -( !2.65p the explicit form of the E-mode polarization can be written: 



AE(ri, r) = -\il- /i2)A'4 - 4/iA'Q - 2Aq ^ 



dlAc 



a^A'u - -^5^Au 



In the case of the magnetized adiabatic mode Au and Aq do not have azimuthal dependence. 
Furthermore a B-mode polarization is only generated thanks to Faraday mixing which has 
been investigated analytically elsewhere [601 EI] and will not be repeated here. The total 
polarization degree coincides with the contribution of Aq, i.e. AQ(rT,, r) = Ap(n, r). It 
follows from Eq. (12.721) that, in Fourier space, the E-mode polarization is 

AE(fc, /i, r) = -dliil - fi')Ap{k, /i, r)], (2.73) 

where Ap obeys 

A'p + {ikfi + e')Ap = Sp{k, fi, r), Sp{k, /x, r) = ^(1 - fi')Sp{k, r). (2.74) 

Using line of sight integration the formal solution of Eqs. (I2.66P and (I2.74p can be formally 
written as [731I7I]: 

Ai(k, /i. To) = r e*^''("-"o)e-^("'"«) ^i(A;, /x, r) dr, (2.75) 
Jo 

Ap{k, fi, To) = r e*'=^("-"°)e-^("'"°) ^p(A;, /i, r) dr. (2.76) 
Jo 

The solution expressed by Eqs. f l2.75p and f l2.76p assumes, implicitly, that the source terms 
can be independently evaluated either numerically or analytically. The approximation of 
tight Coulomb coupling will now be consistently used. The large-scale magnetic fields will 
then affect electrons and protons whose evolution can be determined in the appropriate one- 
fluid limit. It should be stressed that the approach discussed here is very similar, in spirit, to 
the various semi-analytic techniques which have been employed (in the absence of large-scale 
magnetic fields) by various authors [75l [76l [77] starting with the pioneering work of Peebles 
and Yu [78]. 



3 Temperature and polarization anisotropies 

Assuming tight coupling between photons, electrons and baryons, the evolution of the 
monopole and of the dipole of the brightness perturbations determine the evolution of the 
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source term in the temperature and polarization ansiotropies (i.e. Eqs. f l2.66l) - fl2.67l) and 
(I2.74P ). The monopole and the dipole obey, in Fourier space, the following pair of equations 

(i/.- A,„)' = fcA„, (3.1) 
l(R, + l)A„r + 2|(iJ, + 1)A„ = ^A,„ + M^|±i)^ + k(ns-4as) ^ ,3 2, 

where Rh has been already introduced in Eq. (11.41) and where fco is the wave-number 
corresponding to diffusive damping, i.e. the wave-number at which diffusive effects start 
being important. To lowest-order in the photon-baryon coupling the diffusive damping is 
simply proportional to the shear viscosity coefficient t] which has been already introduced 
in Eq. (12.381) . More precisely, to lowest order in the tight-coupling, fc^^ = r]/[p^(l + R^)] 
where 1] has been defined right after Eq. (I2.38P and is proportional to the photon mean 
free path. The estimates based on shear viscosity can be improved by going to higher order 
in the tight-coupling expansion and by further refining the estimates depending upon the 
explicit values of the ACDM parameters. In particular, for typical values of the parameters 
close to the best-fit provided by the ACDM model the values of and £d (i-e. the diffusive 
multipole) can be estimated as 



VrR* + 1 - Jr^ 



The (comoving) angular diameter distance at has been rescaled, in Eq. (13.31) as 

DAiz.) = r^rr Mz*)- (3.4) 
V "MO-no 

Furthermore, always in Eq. (13. 3p tr* is the ratio of the radiation and matter energy densities 
at , i.e. 

r-R* = — ^ = — = 4.15x10 ujy^ [—]. 3.5 

where, following the customary notation, uu = hlQuo- The numerical content of Eqs. 
(I3.3p - (l3.5p is fully specified in terms of whose explicit form can be written as 

= 1048[l + (1.24x 10-=^)cJb°-^^^][l+^i^J'], (3.6) 
0.0783 cub°-^^^ _ 0.560 

Equations (I3.6p - (l3.7p imply z^, = 1090.5 in excellent agreement with the estimate of [Il[2|,[3], 
i.e. = 1090.51 ± 0.95. The evolution of the monopole and of the dipole can be determined 
from the WKB solution of Eqs. and dS^D, i.e. 

_ 

Aio(fc, r) + (P{k, t) = C{k, t) + ^Mik, t) cos [A;r,(r)] e ^, (3.8) 



Aii(A;,r) = 4/'A^(A;,r)sin[A;rs(r)]e '^d , (3.9) 



19 



where C{k,T) and Ai{k,T) are fixed once the initial conditions of the Einstein-Boltzmann 
hierarchy are specified. In what follows, as already mentioned, the initial conditions shall 
correspond to the magnetized adiabatic mode. In Eqs. (13. 8p and (13. 9p rs(r), is the sound 
horizon 



dr 



(3.10) 



/3[i?b(r) + 1] 

whose explicit form will be determined as a function of z^,. The explicit solution of Eqs. 
(12. 9p . (12.101) and (12. lip for the matter radiation transition implies that a{x) = aeq[x^ + 2x] 
with X = t/ti. This also means that: 



+ 1 



'1 



(3.11) 



where ti = 2y (acq/^Mo)/-f^o- By definition of baryon-to-photons ratio (see, e.g., Eq. (11.41) ) 
we have that Rhix) = -Rb*^R*(a;^ + 2x). Thus, defining Eq. y = x + 1, Eq. (I3.10p becomes 
easily 

r"^*/"^! dx 



rs(r*) 



Tl 



V^Rh*rR^ 



Rh*rR*{x'^ + 2x) + 1 

{r./ri) + l dy 



yo 



'1 - Rh*rR^ 

Rh*1"R* 



(3.12) 



which can be integrated via a further change of variables (i.e. y = yosinhw); the result is: 



V3R 



b* ^R* 



arcsinh 



in/n) + 1 



yo 



1 

yo 



(3.13) 



Since, by definition, arcsinh(a) = In [a + + 1] and ti is given after Eq. (13. lip the sound 
horizon at t^, is given by: 

2 1 



1 



HoV^J3Rk,{z, + 1) 



In 



yi + -Rb* + y/Rh^y/l + Tr^ 
1 + VrR* Rh* 



(3.14) 



Having determined the monopole and the dipole by solving Eqs. (13. ip and (13. 2p the polariza- 
tion observables depend chiefly upon the value of the dipole as it arises to lowest order in the 
tight- coupling expansion. However, as it was already observed long ago [73J, the flrst-order 
tight-coupling estimate is not satisfactory from the numerical point of view and must be 
improved. Following this logic, Eqs. ( 12.66P and ( ]2.74p imply the following relations between 
the multipoles of the intensity and polarization brightness perturbations: 



A^o - -[Ap2 + Ai2 - Apo] = -kApi, 



12 



a; 



P2 



9 1 

— Ai2 - — (Apo + Ap2 
L^Ap2-^(Apo + A12) 



-kA 



13 



-kA 



II ! 



3 2 
--/cAps + -/cApi. 



(3.15) 
(3.16) 
(3.17) 
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Summing up Eqs. (13.151] . fl3.16p and (13.171) and recalling that, by definition, Sp = (A12 + 
Apo + Ap2), Eqs. ([3I5D-([3I7D imply 



10 



-e'Sp = k 



-An 



A 



PI 



Ap3 + A 
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(3.18) 



The result of the solution of Eq. (I3.18P turns out to be more accurate than the lowest order 
tight-coupling result. Indeed, neglecting Api, Ap3 and A13 (which are all smaller than An) 
Eq. (13.181) can be formally integrated: 



,T0)/10 



Sp{k,T) = -/te3^(-'-o)/io r dT'An{k,T')e-^<^'^ 
5 Jo 

which also implies that 

Ap{k,fi,To) = -0.515 (A; - fi^)e"'^'^^'-^°^An{k,n) 



(3.19) 



(3.20) 



The coefficients and can be determined in terms of Eqs. (I3.8p - (l3.9p and (13.20p 
following the standard techniques. More specifically the coefficient a\J^ can be expressed as: 

1.2 



a 



(T) 



(27r)3/2 



(Aio + (p)jiix) + 3 An 



dje 
dx 



(3.21) 



where a; = fc(ro — r*) and where it{x) are the spherical Bessel functions |6Hl |69] of argument 
X. In Eq. (I3.2ip kt = v^/c* arises from the integration over r of the Gaussian visibility 
function. The coefficient turns out to be: 



3 



(£-2)! 



47rv/2£Tl j d^kx^ [{l + diy]je{x) J^' IC{T)Sp{k,T)dT. (3.22) 



4(27r)3/2^ {i + 2) 

In Eqs. (13.211) and (13.221) the following two results have been repeatedly used: 



dn Y,* (n) e 



2\2 —ifixi 



^Ani-tYVWTljfXx), 
dl[il + dire-^n = -a + dl)x'e-^>^^ 



(3.23) 
(3.24) 



Furthermore, in Eq. (I3.22p . the equation of the spherical Bessel functions [68l 169] has been 
repeatedly used. Notice that, in Eq. (13.220 the integral over r of the visibility function can 
be simplified by using, for Sp{k, n, r), the expression of Eq. (13.190 and by performing exactly 
the same integral leading to Eq. (13.200 . 

In the present discussion we are interested in the scaling properties of the correlation 
functions over relatively small scales where simplifying expressions for the Bessel functions 
can be used. In this limit, the reionization effects can be parametrized as follows. In Eqs. 
(I2.75P and (12.761) the integral over r can be separated in two distinct contributions. For sake 
of concreteness consider the polarization integral which gives 



Ap(A;,/i,To) 



TO 



rfr/C(r)e-*^^'^p(A;, /x, r) + / rfr/C(r)e~*'^^^p(fc, /z, r). 
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(3.25) 



The first term at the right hand side Eq. fl3.25p is the most relevant for i ^ 20 and it is 
given by e~'''^'=Ap(A;, /i, tq) where Ap(/c,/i,ro) is the value of the polarization in the absence 
of reionization. For small i the second term in Eq. fl3.25p is the most relevant [67] and it 
leads to supplementary peaks in the angular power spectra (i.e. the so-called reionization 
peaks). Within the approximations of this section, the integrand of the second term in Eq. 
f l3.25p is simply proportional to the quadrupole of the intensity which can be evaluated, for 
k < /cd, as [67] 

Ai2(A;,/x,rre) = {C{k,n) + ^/c~^M{k,T^) cos[krs{n)]}j2ixre) (3.26) 

where j2(a^re) is the spherical Bessel function for £ = 2 and where = [k{Tre — t*)]. The 
reionization peaks arise, roughly, at the first peak of j2{xy-e), i.e. for — 2. 



4 The basic integrals 



4.1 Generalities 

The considerations of the previous section depend upon two sorts of scales i.e. 

• damping scales, (e. g. 4, Cre---) which control the falloff of the temperature and 
polarization angular power spectra; 

• oscillatory scales (e.g. Ea) which control the structure of the peak and depths in the 
TT, EE and TE correlations. 



The thermal diffusivity multipole £d, already introduced in Eq. (13. 3p . can be estimated using 
the best fit to the WMAP Syr data alone [HEIE]; Eq. ([S3]) leads to = 1422.08. The finite 
thickness of the visibility function leads to an effective multipole which can be estimated as 
A = V^/ikQa^). Again using the estimated thickness in the visibility function the WMAP 
Syr data allow to estimate £t = 1211.22. The typical scales £■£, and it can be combined in 
what is often called Silk damping scale, i.e. 



(?2 + ^2 ' 



1 



(4.1) 



In the case of the numerical values listed above is = 922.09. The oscillatory patterns in the 
angular power spectra are determined by the acoustic multipole, i.e. 



27r(iA(2* 



A/3i?b* TT A/i~FT(iA(2* 



In 



Vl+fib* + A/(l+rR,)fib, 



~ U03 



In 



i+Vm.-Rb. 



(4.2) 
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The first equality of Eq. (14. 2p is just tlie definition of tlie acoustic multipole wfiile tfie second 
equality uses a more explicit form of the (comoving) angular diameter distance. Note that, 
in Eq. (14.21) . the (reduced) angular diameter distance dxiz^:) goes asymptotically to 0.89 
for z^, > 500 and for the standard values of the cosmological parameters. In the explicit 
expressions of some integrand it will prove useful to have an explicit expression also for 
7a = I^K- indeed, in various oscillating factors, the combination arises naturally. Thus 
Eq. (14. 2 p also implies 



7a = ^ = , = In 

rfA(2*)V3i?b*(^* + 1) 



yTT:R^ + ^(l + rR,)i?b 

1 + aArJ^ 



(4.3) 



According to Eq. (14. 2p . £a = 301.57 while the WMAP 5yr data imply that the acoustic scale 
at decouphng is given by £a = 302. 08*^0^84. The acoustic multipole can be compared with 
the equality multipole, i.e. 



.eq - v^SfiMo^^oy^^cq + \B p^z^^ = 2\pl .Jz^^^Ji^ = 2^2 ^ ^ d p,{z^) , (4.4) 

where we recalled that Tieq = ^eq-f^eq = V 2VLy[o (oo / c^eq) , Combining the last equality of 
Eq. (14. 4p with Eq. (13. 5p we get £eq = 439.057 Y/cJM(iA(-2*) which also implies, in the case 
of the WMAP 5yr parameters, £eq — 136.95. In terms of £a the position of the first three 
Doppler peaks can be obtained approximately obtained from = ^a{.^ — V^m) where 

£i = £a(1-¥^i), y^i = 0.267 (^y'', (4.5) 



4 = ^a(2-^2), ^2 = 0.24 (^^j , (4.6) 

4 = ^a(3-¥^3), = 0.35 (4.7) 

The values of ii, £2 and £3 are deduced in the case Ug = I. When ng 7^ 1 the positions are 
shifted as im ^ (^m + 

A£i = 0.13|ns-l|£i, A^2 = 0.33 Iris -11^2, Afs = 0.61 |ns - 1I4- (4.8) 
In the vanilla ACDM and for the WMAP 5yr best fit we have that 

£1 = 219, £2 = 535, h = 814. (4.9) 
which is approximately what could be obtained from Eqs. (I4.7p - (l4.8p 

4.2 The angular power spectra 

Using Eqs. (I3.2ip and (13.220 into Eqs. (I2.60p . (I2.6ip and (12.620 a more explicit expression 
of the angular power spectra can be obtained: 

Cr'=4./f|l|Ar'(*,r„)|^ (4.10) 
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where the following quantities have been introduced: 
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lA 



(TE), 



(4.11) 
(4.12) 

(4.13) 
(4.14) 
(4.15) 



It is practical to adopt the following general parametrization for the three relevant power 
spectra of the problem 

\ ris-l / \ 2(nB-l) 



Vnik) = A 



n 



Vnik)=£^ 



VM = rBVnik), 



(4.16) 



where An denotes the amplitude of the curvature perturbations at the pivot scale k^, E-q 
denotes the amplitude of the power spectrum of f^B (see also [IHl HH [60j; re denotes the 
ratic@ between the power spectrum of ub and the power spectrum of fie at the same magnetic 
pivot scale k\^. To leading order and tb are independent upon the wave-number. There 
are however corrections which imply that £^b and tb do depend upon the wave-number. If 
£^(/c) and r^ik:) the form of the integrals listed below as well as the related discussion slightly 
changes but the explicit results are more cumbersome and will not be reported here. 

Using Eqs. ([31]) and ([31]) into Eqs. fHl3]l . fITO]) and fITO]) the exphcit form of the 
temperature and polarization observables can be derived. Since some of the subsequent 
expressions are rather lengthy, the following rescaled amplitudes will be defined: 
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(4.17) 



where kp = 0.002 Mpc ^ is the pivot scale of curvature perturbations and, as already men- 
tioned in section [2], /cl = 1 Mpc~^ is the magnetic pivot scale. 



4.3 Temperature autocorrelations 

The temperature autocorrelations are hereby written in terms of four basic integrals, i.e. 




= X''>{i,Q+l''\e,is)+Vs)\£Js)+V^''{i,iuis). (4.18) 



^It is often practical to assign ratios of power spectra at the same pivot scale; this is what happens also 
when assigning tensor power spectra in standard CMB studies. 
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Each of the terms appearing in Eq. fl4.18p contains three contributions proportional, respec- 
tively, to Q-R-R, Qbb and to Q-rb whose explicit form can be written as: 
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(4.19) 
(4.20) 
(4.21) 

VS(^, ^s, A) + V^1^(£, 4, it) + cos/5(vS + V^'i)(£, is, it), (4.22) 

where cos f3 parametrizes the correlation between the purely adiabatic and the purely mag- 
netized component^. The terms appearing in Eqs. fl4.19p are expressible as 
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bb(^? it 



vS(AA) 



'■ Q-R-R 

QBB^'^""-'^/S(^,A,2nB-l), 
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The basic integral appearing in Eqs. fl4.23l) - fl4.25l) is given by 0: 



/S(^,4,n) 



n— 3 



--Lx{w,i)LY{w,- 



dw. 



(4.23) 
(4.24) 

(4.25) 
(4.26) 



In Eqs. fl4.26p the functions Lxiw, i) and Ly{w, i) account for the contribution of large-scale 
magnetic fields to the tight coupling solutions and also depend upon the resolution of the 
calculation, i.e. upon i^^ax (which denotes the maximal multipole at which the calculation 
is trustable). The four functions which enter Eq. (14.261) as well as the other seven integrals 
which will be discussed below are: 



a-R - P-R^n {wqi), 
--an + Pn In (w qe) 



Lb{w, £) = aB- 13b In {w qi), 
Mb{w, i) = aB + Pb In {w qe) 
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(4.27) 
(4.28) 

(4.29) 
(4.30) 
(4.31) 
(4.32) 



^This correlation arises, in explicit models, because magnetic fields are produced during some stages of 
inflationary expansion [23] (see also [24l [25l [26l [27] ) . In a model-independent perspective the correlation 
between different components should also be considered in full analogy with what happens for entropic initial 
conditions [SD] [H] [51 [H [SI] . 

^"^It is relevant to point out that the arguments of the integrals contain the multipole, the diffusion scales 
and the relevant spectral index. These are the basic quantities which define the eight basic integrals which 
will now be listed. 
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where i?b is the baryon to photon ratio at the recombination and qi is given by: 



^^-(2oorfA(..))\/fTI- ^^-^^^ 

Concerning the notations of Eqs. ( I4.29p - (]4.32p we remind, as defined after Eq. (14.160 that 
'"B = Vn/Va- is, therefore, the ratio between the power spectrum associated with Qb and 
the power spectrum associated with ctb. 

Since the aim of the present analysis is to have analytic estimates of the modifications 
induced by large-scale magnetic fields especially at small angular scales (i.e. in the limit 
£ 3> 1). In the latter limit the spherical Bessel functions je{x) can be approximated in their 
large-order limit and the acoustic multipole fixes the oscillatory structure of the angular 
power spectra: 

iii + ^ iii + lf-^^)§^ ^ I .\ , x = wl, (4.34) 



where /3(x, £) = y/x'^ ~ P — i arccos {i/x) — j [68l[69]- Recall that, often, changes of variables 
are required to evaluate the integrals. In particular, a practical choice is: 



o , , dw 2dy , ^ 

w^y' + l, dw^2ydy, ^ 4.35 

Vw^ - 1 Vy + 2 

The change of variable = + 1 is also possible in some cases and it leads to a simpler 
structure of the integrands, in some cases. In spite of the change of variables, the numerical 
values of the various integrals do not change. At the same time, since the integrals will be 
evaluated numerically, the time of evaluation can also change as a function of the algebraic 
form of the various integrands. 

The contribution labeled by I^2)^\£,is) in Eqs. (14.181) and (14.201) leads to the following 
explicit results 

y^k^Js) = Q7^7^^-^4'4(£,4,n.), (4.36) 
Vi(^,4) = Qbb e^^^-'^lSk^,is,2nB - 1), (4.37) 
n„„^^^^i^^r(2)/« « + 2nB - 1 



mi^Jt) = QnBi^^I^i,[iJs, / j, (4.38) 

where the second basic integral appearing in Eqs. (I4.36p - (l4.38p can be written as: 

I^xk^,iun) = ^J^^W4w,c,^)w''-'^MxiwJ)My{wJ)e~"^k^'"'dw. (4.39) 

For practical convenience, the two functions W±{w,Csh) are introduced, respectively, in Eq. 
fO^ and in Eq. flCTj) : 

W,Kc..) = ^t<ii||l£i^, (4.40) 

V^f — 1 



26 



where Csb is the photon-baryon sound speed aheady introduced, for instance, in Eq. fll.4p . 
The third contribution appearing in Eq. fl4.18p . i.e. X\ is determined by the terms 

appearing in Eq. (14.211) whose exphcit expressions are: 

V4t(^,^s) = Qnnr^-'lgk^,is,n,), (4.41) 
^Sk^Js) = QBB^'("^^''4'iJ(^,4,2nB - 1), (4.42) 
;(3)^. n__^^^^i^r(3)^« « + 2nB - 1 



V^^^(^, it) = QnB ^^^4b [i, 4, j . (4.43) 

Recalhng Eq. (14.401) the basic integral appearing in Eqs. (I4.41I) - (I4.43I) is given by: 

jjf|.(^,4,^) = 2/ W-{w,c,^)w''-^ cos{2^Aiw)MxiwJ)MY{wJ)e ^ dw, (AAA) 



where 7a has been introduced in Eq. (14. 3 p and >V_(w,Csb) is defined in Eq. (14. 40 p . Finally, 
the fourth basic term of Eqs. (I4.18P and (I4.22p is completely specified by the four expressions: 

V4t(^,^s) = Q7^7^^-^J4'4(£,£s,A,r^s), (4.45) 
VSJ(^,4) = QBB^'^"«~'^4lJ(^,4,A,2nB - 1), (4.46) 

VS(£,£0 = 6.B^"^4t(^,4, "- + ^;^"7 (4.47) 
The basic integral appearing in Eqs. ( 14.4ip -( l4.43p is given by: 

I^fiyii,isJun) = 2 ^L^==cos{^dw)Lx{wA)My{wA)e'^s dw. (4.48) 

Ji \Jw^ - 1 

Equations (14.261) . (14.391) . (I4.44p and (I4.48P define the primary form of the integrals deter- 
mining the temperature autocorrelations. In what follows the EE and TE correlations will 
be more specifically studied. 

4.4 E-mode autocorrelations 

Within the same logical scheme already employed in the case of the TT correlations, the EE 
angular power spectra of Eq. (14.110 can be written in terms of two (further) basic integrals, 
i.e. 

^(EE) ^ ^(EE) _ ^(EE) _ ^4_4g^ 

Both the EE and the TE angular power spectra are suppressed with respect to the TT 
correlations. It is therefore useful to define the quantity 

Ar(^E)(^^ ^ Q_^32 ^^^^j2 ^ ^)2(^ - 1)(£ + 2) r^ (4.50) 

which is independent of I in the range of multipoles where the calculation can be trusted 
(i.e., in practice, i > 20). The explicit form of the integrals appearing in Eq. ( 14.491 can be 
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written, in full analogy with Eqs. (K1M - (KT3\ . as 



where 



4T^(^,^d) = vSU^Jd) + vSiiiJn) + 2cos/3vS(£,£d), (4.51) 
^le^\^jB) = vS(£,£d) + vS1(^,^d) + 2cos/3vS(£,£d). (4.52) 



V4'i(^, iu) = QnnAf^''^\(^, c^I^U^, £d, n,), (4.53) 

vS(^,^d) = QBBAr('^''n^,^*)^'""d4'i(^,^D,2nB - 1) (4.54) 

vS(£,£o) = Qn.Af^--^ a, a.) ff, + " ^ (4.55) 



where 

lfy{l,i^,n)= ^-^ Mx{w,i)MY{w,i)e ^ dw. (4.56) 

The three terms defining x|™^(£, ^d) are: 

vSk^, in) = QnnM^^^\i, a^)t^-' c^jSizi^^, £d, n^), (4.57) 
vSkiJu) = QBBAr(E^H^,^*)^'"""'d4'i(^,^D,2nB - 1), (4.58) 
V^lUM = Q.BAr(-)(£,a.)£-^4/S(£,£.,!^i±f^). (4.59) 

The sixth basic integral appearing in Eqs. fl4.57l) - fl4.59l) is 

/(«)«, I,, n) = r "^P^f'';-' M,(™, l)My(r.. Oi^i'"'^ (4.60) 
Ji \/w^ - 1 

Equations fl4.56l) and fl4.6UI) represent the primary form of the two basic integrals determining 
the polarization autocorrelations. From the purely algebraic point of view the EE angular 
power spectra have a single periodicity is insofar as Eq. (14.561) has an integrand which 
is not oscillating while the integrand of Eq. fl4.60p depends on a single oscillating term. 
Conversely, the TT angular power spectra are given by the weighted superposition of the 
integrals appearing in Eqs. fl4.26p . (14.391) . (I4.44p and (14.480 whose corresponding integrals 
do not depend upon the same oscillating term. The single periodicity of the EE angular 
power spectra will, have, in the present context, interesting consequences. 

4.5 Temperature-polarization cross-correlations 

The last angular power spectrum considered here is the one arising from the cross-correlations 
between temperature and polarization, i.e. the TE power spectrum leading to the following 
integrals 
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where 



In full analogy with what has been done in the case of the EE correlations (see Eq. fl4.50l) ) 
it is practical to define 

Ar(™)(£^^^) = 0.515 fcoa, (^ + 1) v'(£ + 1) (£-!)(£ + 2) r 5/2. (4.62) 

Consequently, the explicit form of the integrals appearing in Eq. fl4.6ip is 

X}™)(£,£s) = V^k^Js) + vSkeJs) + cos/?[vS(£,£s) + vSk^Js)], (4.63) 
X(™)(£, is, ^d) = VS(£, is, in) + vS^ii, is, in) + 2 cos/3vS(£, is, in). (4.64) 

vSkii, 4) = Q7^7^Ar(™)(£, a.)r= cfll^U^, is, n,), (4.65) 
vSk^, 4) = QbbAT^t^H^, a.)^'"«-^ 4, 2nB - 1) (4.66) 

vSUi,is) = Q.BAr(-)(£,a.)cf £-^/S(£,£s,^^^^^^). (4.67) 

The integral appearing in Eqs. fl4.65l) - fl4.67l) 

./i - 1 

The last bunch of terms contributing to Ig^^\i,is,iD) is given by 

V4t(^, (s, iB) = QnnJ^^^''\i, a*)r= 44']e(^, 4, ^d, n.), (4.69) 
V^?(^, 4) = QBBM^^'^Ki, fT*)^'"^"^ d/S(£, 4, ^D, 2nB - 1), (4.70) 

vS(£,4,^d) = Q7^BA^(T^^^,^*K"^d4'i(^,4,^D^^^^^^), (4.71) 

where 

i^kej,M = i r!i!i(»£!:!rf..M,(u,,QA/,(..,Qe-'t"f (4.72) 

2 Ji — 1 

The integrals of Eqs. (14.681) and fl4.72p give the last pair of primary integrals. The results 
obtained in the present section allow for an explicit evaluation of the TT, EE and TE 
angular power spectra. The following section is devoted to the derivation of a number of 
scaling relations which are the magnetized analog of the standard scaling relations which 
constitute the basis of any sound strategy of parameter estimation. 



5 Scaling properties and form factors 

The 8 basic integrals derived in section H] can be exploited to study the deviations induced by 
the ambient magnetic field on the CMB observables. In the present section the semi-analytic 
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Figure 3: The semi-analytic results for the polarization autocorrelations are illustrated in 
the absence of magnetic fields. In the plot at the left, the explicit result for the integrals 
of Eq. fl5.5p is reported. In the plot at the right the explicit result of Eqs. fl5.6p - fl5.7p is 
confronted to the WMAP Syr best fit. 



results will be confronted with the numerical estimates. The purpose will not be to touch 
upon all the possible themes of the analysis but rather to mention only some of the most 
notable aspects which emerged from an exhaustive study of these matters. The polarization 
autocorrelations are sensitive to 2 out of 8 basic integrals and, as previously discussed (see, 
e.g. discussion after Eq. fl4.60p ). they have the simpler periodicity. They also depend upon 
£d since the integral over the optical depth allows for an explicit integration of the source 
term (see Eqs. f l3.15p - fl3.19p ). The present section is organized as follows: in subsection 
15.11 the EE angular power spectra will be discussed and the semi-analytical results will be 
compared with the numerical evaluation. In subsection 15.21 the semi-analytical results for 
the TT and TE correlations will be illustrated. Finally, subsection 15.31 will be focussed on 
the scaling properties of the temperature and polarization autocorrelations. 



5.1 EE angular power spectra 

In the absence of any ambient magnetic field, Eqs. fl4.49p and fl4.5ip - fl4.52p lead to the 
complete expression of the EE correlation which can be written as 

^(EE) ^ jg^)(£,£D)-X(7\£,^D), 
X(TH^,^d) = V4^4(£,£d), X(f)^H£,£D)=V4^4(£,£D). (5.1) 

Bearing in mind the explicit form of the different contributions, Eq. fl5.ip becomes 

^(EE) ^ -^(EE) _ ^^^^ ^ ^ 2)r-3[/S(£, £d, ns) - 4'i(£, £d, ris)], (5.2) 
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a!^''''^ = (0.132) {koa^r (^)"" cle-'^-AnT^, (5.3) 

where A is the rescaled amphtude grouping all the factors which are independent on the 
multipole i. If Us = 0.963 and ere = 0.087 (as in the Syr best fit to the WMAP data alone), 
Eq. ([E3D implies that A^^^^ = 4.25 x 10"^ (/iK)^. The factor (fcocr,) can be estimated 
within the WMAP data and it is given by 1.43 x 10~^. The latter figure arises, as discussed 
after Eq. fl2.35p . by computing /S.t^/tq where Ar* is the thickness of the last scattering 
surface in conformal time. The integrals of Eq. (15.21) appeared already in Eqs. (14.561) and 



n„ =1 .5, B = 1 nG H = 71 .9, £2^ „ =0.0441 , £2 =0.21 4, n. =0.742, e = 0.087 
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Figure 4: The semi-analytic results for the polarization autocorrelations are illustrated in 
the presence of magnetic fields. In the left plot the full line denotes the numerical result 
while the dashed line denotes the analytic approximation. In the right plot the magnetized 
result is compared with the WMAP best fit. In this and in the following plots /? = 0. 

(I4.60p and their explicit expressions, for the case at hand, is: 

4t(^,^D,n.) = ^rf«;^^M^(^,£)e-^(^/^-)^-^ 
Ji \Jw^ — 1 

/S(^,^D,ris) = / dw cos M^w, i) e~^^"'-^'^\ (5.4) 

Ji — 1 

It impossible to change integration variable in Eqs. (15. 4p . By positing w = + 1 we do 



f EE) 

^^It should be noticed that the expression for A is dimension- full since the result has been multiplied, 
as customary, by T^q = (2.725 x 10^)^ (mK)^ where T~^o is the inferred value of the CMB black-body spectrum 
in units of /iK. 

^^The change of variable w'^ = -f 1 is also possible and, in this particular case, will lead, of course, to the 
same results. In the case of other integrals, however, mathematically equivalent change of variables might 
lead to different evaluation times of the corresponding numerical integrals. 
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4%(^,^D,ns) = dy ^y cos[27A^(y^ + l)]M^(y,£)e" 'V \(5.5) 

The integrals of Eq. flS.Sp converge rapidly and can be estimated, for instance, with numerical 
techniques; the final result can be expressed in a closed form for £ > £i as 



Gf ^) = + - 6e cos [27a(^ + ^i)]} e-2(^/^-)^ (5.6) 

:4^''^^ = 4.476 X 10-^ (0.0354)--^ ( 2 4lt'\o-9 ^^''^^ 
aE = 1.67, 6e = 3.38, h = 65, (5.8) 



where £1 appears because the analytic derivations of the previous sections assume a large- 
order expression for the spherical Bessel functions. Concerning Eq. (15. 7p few comments are 
in order: 

• Eq. (15.71) assumes the simplified treatment of reionization which has been spelled out 
in Eqs. ( I3.25p - (l3.26p and which is less accurate for low multipoles (i.e. in the region 
of the reionization peaks) than for large multipoles; 

• oe and 6e are, respectively, the form factors coming from the integral X^]j(£, f^, '^s) 
and from I^]j(£, Id, Us); 

• £d (i-e. the diffusion damping scale) is given by Eq. ( 13. 3p and has been also discussed 
prior to Eq. (14.11) in connection with the estimate of Silk damping; 

• the numerical value of vA*-™-* follows from the pivotal value of kp (i.e. 0.002 Mpc^^) 
and by computing ko from the (comoving) angular diameter distance of Eq. (13.41) : 

In terms of the values of the cosmological parameters obtainable from the WMAP 5yr best 

fit [DEI E] 

{um, Wc, cJb,cJA, ho, n„ ere) = (0.1326, 0.1099, 0.02273, 0.385, 0.719, 0.963, 0.087), (5.9) 
the values of the derived parameters of Eqs. (I3.3p - (l3.6p ar 

[z^, Csb(2;*), /^a(^*), ^a, £d, a, 4] = 

[1099.5, 0.451, 14110.8 Mpc, 301.578, 1422, 1211, 922]. (5.10) 

The results for the polarization autocorrelations are illustrated in Fig. [31 In the plot at 
the left I^]j(£, £d, ""-s) and I^]j(£, £d, "^s) are reported, respectively, with the full and with 
the dashed fines. In tfie plot at tfie rigfit of Fig. [3] tfie analytic result for G^™^ (dasfied 



^•^Different best- fit parameters, obtained by combining CMB data with other data sets (e.g. [79l [80] ) 
lead to different values of the derived parameters which can be however computed always using the general 
formulae of the previous sections. 
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line) is compared with the WMAP Syr best fit (full line) holding for exactly the same set 
of parameters (i.e. Eq. (15. 9p ). In Fig. H] the analytic results for the magnetized polar- 
ization autocorrelations are compared with the numerical results. In both plots of Fig. H] 
the numerical result is reported with the full line. The dashed line denotes the analytical 
approximation (plot at the left). The dot-dashed line denotes the WMAP Syr best fit (for 
the same value of cosmological parameters). In both plots the correlation angle has been 
chosen as /5 = 0. In summary we can therefore say that Figs. [3] and H] show that, in spite of 
the different approximations, the analytic result is in fair agreement with the numerical one. 
Finally, the numerical results illustrated in Figs. H] and [5] follow from an improved version of 
the approach already mentioned in the introduction [481 149| which is based on a modification 
of CMBFAST tSUES] (which is, in turn, a modified version of Cosmics pjlM]). 

5.2 TT and TE angular power spectra 

The TT and TE correlations share similar features from the point of view of the analytic 
results discussed here. The periodicities of the TT and TE angular power spectra arise as the 
weighted interference of the periodicities of the monopole and of the dipole of the radiation 
field. The TT correlations, have been partially discussed with a similar semi-analytic method 
in IIH] (first paper) and corroborated by subsequent numerical estimates (second paper of 
[H]). The improved analytical understanding developed in the present paper allows for a 
better assessment of the accuracy of the results. In Fig. [5] the results for the TT angular 
power spectra are illustrated. In both plots the full lines denote the analytical estimate while 
the dashed lines represent the numerical result. In Fig. \5[ from left to right, the magnetic 
field intensity and the spectra index increase. In both plots of Fig. [5] the dot-dashed lines 
denote the WMAP Syr best fit. The results illustrated in Fig. [5] are representative of a 
general trend which has been observed also in other cases and can be summarized as follows: 

• the analytic result for the TT correlations stemming from the basic integrals studied 
in this paper always underestimates the numerical result; 

• the analytic result becomes progressively inaccurate as the field strength increases 
above 10 nG. 

Various other examples show that the polarization autocorrelations (i.e. the EE angular 
power spectra) are better captured by the analytical results, as already shown in Figs. [3] 
and m The temperature-polarization cross-correlations (i.e. the TE angular power spectra) 
share the same levels of accuracy of the TT correlations and they are illustrated in Fig. O 
As in Fig. O the dot-dashed lines denote the best fit to the WMAP Syr data alone. For 
-Bl < 10 n G the analytic results are rather accurate for the TT, TE and EE angular power 
spectra. In the case i?L > 10 nG the results become progressively less accurate in the case 
of the TT and TE correlations but remain reasonable for the EE autocorrelations. 
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5.3 Envelopes and wiggles 

Having tested the accuracy of the analytical results, the handiness of the approach developed 
in the present investigation resides in the determination of the scaling properties of the var- 
ious correlation functions. To proceed in this direction, the idea is to compute numerically 



H = 71 .9, = 0.0441 , £1 =0.214, £2. 0.742,6 = 0.087 




H„ = 71 .9, a „ = 0.0441 , £2 „ =0.21 4, £2. =0.742, E =0.087 




Figure 6: The semi-analytic results for the temperature-polarization cross-correlations are 
illustrated always in the case f3 = 0. 



the ratios 

(TT) _ Gf^H^B,^L) p(EE) _ Gf^\nB,Bi^) 

- (TT) ' ~ ^(EE) ' l^'-L-LJ 

where, by definition, G^^'^'^ and G^™'' denote the angular power spectra in the absence of 
ambient magnetic field. The same procedure can be carried on also for the TE correlations. 
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Figure 7: The magnetized form factor from the TT correlations for different values of the 
magnetic field background. The correlation angle is fixed to /3 = 0. 



However, since the TE corrrelations are not positive definite and pass through zero, the 
resulting plots are not as revealing as the ones obtainable from the TT and EE angular 
power spectra. Figures [7] and [8] illustrate, respectively, R} and R} for different values 
of the magnetic field intensities and for different values of the spectral indices. The ratios 
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Figure 8: The magnetized form factor from the EE correlations for different values of the 
magnetic field background. The correlation angle is fixed to (3 = 0. 

R^'^'^^ and R^^^^ represent an effective numerical diagnostic of the possible influence of a 
putative magnetic field. Indeed the diffusive scales, the thickness of the last scattering 
surface the optical depth at reionization are only mildly sensitive to the presence of an 
ambient magnetic field. The analytic structure of the angular power spectra computed in 
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section H] suggests the possibility of factorizing the effects of a putative magnetic field into an 
appropriate form factor which is, incidentally, illustrated in Figs. [7] and [8] with a dashed line. 

(XY) 

On a qualitative ground the form factor represents, in some sense, the average 

of i?r • On 

a more quantitative ground the results of section H] suggest that -R^^^'* can be factorized into 
the product of a non-oscillating factor (i.e. F^(nB,-BL)) and of an oscillating contribution 
(i.e. Of^^\ns, ne, Ati,Bi^)). Therefore, within the notations followed in the present paper, 

(XY) 

we will have that R\ can be written as 

R\'''^\nB,Bi^,n,,An) = Fi{nB,Bi^)Of''\n,, ub, An, ^l). (5.12) 
where X,Y = T, E. The result for the magnetic form factor is 

Fe{nB,B^) = l + ai(^^y Ji{n,,nB,i) + a2(^^yj2{n,,nB,f3,i). (5.13) 

Within the set of parameters given by Eq. (15.91) the constants ai and a2 are given by 

ai = 1.393 X 10^^ as = 1.952 x 10-^ (5.14) 
while the functions j7i(ns,nB,^) and 1/2(^8, ne, /5, ^) can be expressed as 

J,{n,,nB,i)={^) (^^j [^) Si(n„nB), (5.15) 



Sn-R — ris — 1 



j2(ns,nB,Af) = cos/?(^^j (^^j (^-j S2(ns,nB). (5.16) 

In Eqs. (15.151) and (15.161) the functions Si(ns,nB) and S2(ns,nB) encode a milder depen- 
dence upon the spectral indices and they can be determined by matching the Sachs- Wolfe 
expression (valid for i < ii) with the results of the explicit numerical integration of the basic 
integrals (vahd for i > ii). The form of Si(ns,nB) and E2(ns,nB) also depend upon the 
regularization scheme of the magnetic energy density and here the explicit expressions will 
be given in the case of blue magnetic spectral indices: 

SiK,nB) = p., ^'^^^~2^.fo ^ - 1). (5.17) 

r(3 - ns)r^(5/2 - ub) 



^2{ns, ub) = 7 ^ yfi^B - 1), (5.18) 



^^''^"3x(3-2x)r2(i/2)' ^^'^^^ 

(5.20) 

The functions appearing in Eqs. fl5.17p - fl5.18p can also be estimated (just in the limit of 
large i) from the analytic expressions of the different integrals, as illustrated above in this 
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section. In the latter case the resulting expression will still have the correct scaling properties 
but the overall normalization will have to be adjusted. Conversely, the advantage of Eqs. 
fl5.15p - fl5.16p and fl5.17p - fl5.18p is that they are immediately comparable to the numerical 
calculation also for small i. In Eqs. (15.150 and (15.160 we have that is = 1- If the integrals 
would just be estimated from their small-scale approximation the putative value of £b would 
be larger and will fix the limits of applicability of the formula. 

As already mentioned, in Figs. [6] and [7] the dashed curves illustrate the magnetic form 
factor of Eq. (15.130 for the different values of the parameters appearing in each plot. The 
structure of Figs. [7] and M is as follows: 

• in the three plots at the left the magnetic spectral index is fixed to = 1.5 while in 
the three plots at the right the magnetic spectral index is fixed to ub = 2; 

• from top to bottom, as indicated in each plot, the values of the magnetic field intensity 
augments from 1 nG to 20 nG. 



According to Figs. [7]and[8l the magnetic form factor of Eq. 15.131 reproduces quite faithfully 
the average of the numerical data and this is what scaling relations can provide, in this 
context. As explicitly shown by Figs. [7] and [8] the very same form factor works both for 
the TT and for the EE angular power spectra. The latter observation demonstrates that 
the factorization of Eq. fl5.12p is not only analytically plausible but it is also numerically 
justified. The last observation brings up a further question: can we also understand semi- 
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Figure 9: The analytic form of the wiggles is compared with the numerical results. 

analytically the structure of the wiggles of Figs. [7] and [HP While in Fig. [9] the wiggles 
exhibit a double periodicity (i.e. a hump is followed by a peak), in Fig. [H] there is a single 
periodicity (i.e. a single peak is followed by a single peak). The difference in the two 
structure is understandable on the basis of the considerations of sections [3] and H) while 
the TT correlations arise as the interference of the monopole and of the dipole ther EE 
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correlations mainly feel the dipole and are, therefore, a cleaner probe where the analytic 
considerations can be more easily confronted with the numeric results. Going back to the 
parametrization of Eq. fl5.12p the analytic structure of the wiggles can be written as 



0(^^)(nB,i?L) 



qe - bE cos [27a (£ + - 6b] 
oe - bE COS [27a(^ + ^i)] 



9.2 X 10" 



nG 



+ 



5l 
nG 



(5.21) 
(5.22) 



where ii, oe and 6e have been already introduced in Eqs. fl5.6p - fl5.7l) . The rationale behind 
Eqs. f l5.21l) - fl5.22l) is rather simple. In the denominator of Eq. (15.211) there is the analytic 
form of the best fit, while in the numerator the ambient magnetic field introduces a phase 
difference. In Fig. [9] the analytic expressions for the wiggles are compared with the numerical 
results. The remaining offsets are within the accuracy of the analytic approach and improve 
on the pure scaling estimate which lead to the derivation of the envelope. The results 
derived here are relevant for the dedicated strategies of parameter extraction which have 
been suggested in [62] . 



6 Parameter space of magnetized CMB observables 

In the present study the values of the magnetic field intensities and of the corresponding 
spectral indices have been taken to be, in some cases, rather large in the sense that the 
selected values lead to CMB observables which are incompatible with the observed ones. As 
already mentioned the largeness of some of the selected values is evident from the comparison 
of the computed CMB observables with the best fit to the WMAP 5yr data alone (see, e.g. 
Fig. ED 

The choice of dealing with some of these extreme values is, in some sense, dictated by 
the logic followed in the present study: for large values of the magnetic fields, the scaling 
properties of the angular power spectra are more transparent and the distortions enhanced. 
We hope it is clear, from the results of the previous section, that indeed, the distortion 
patterns scale with the amplitude but their morphology remains unchanged. This is, after 
all, closely related to the intuitive notion of scaling (see, e.g. Figs. [7]and[H]). 

At the beginning of Section [2] it has been mentioned that, for instance, the values 
(ne, -Bl) = (2, 10 nG) are excluded, in a frequentist perspective, to 95 % confidence level. 
In what follows the latter statement will be made more quantitative by deriving and by 
discussing the relevant exclusion plots in terms of the spectral index ub and of the mag- 
netic field intensity i?L- Let us, first of all, convince ourselves that, indeed, the model 
{ub, Bi^) = (2, 10 nG) does not correctly reproduce the data. In Fig. [TD]the full line illus- 
trates the best fit to the WMAP 5yr data alone while the dashed line is the numerical result 
for the case (ne, -Bl) = (2, 10 nG), when all the other parameters are fixed as in Eq. (15. 9p . 
If Fig. [To] the data points (with the corresponding error bars) refer to the binned data both 
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Figure 10: The TT and TE angular power spectra are illustrated in the case {n^, 5l) = 
(2, 10 nG) (dashed lines) and for the best fit to the WMAP Syr data alone (full line). In 
both cases [3 = for consistency with the analysis of the previous sections. 



for the TT and TE angular power spectra. The binned data contain 34 (effective) multipoles 
in the TE correlation and 43 (effective)multipoles in the TT spectrum. Following the usual 
habit, to make the plots more readable, the binned data points have been included. Con- 
versely, the unbinned data (which are the ones used in the following analyses) contemplate 
all the multipoles from i = 2 to i = 1000 both for the TT and for the TE (observed) power 
spectra. Finally, we included the TT and the TE angular power spectra since they are the 
best measured spectra in the context of the WMAP Syr data. 

Having established that the parameters {ub, Bi) = (2, 10 nG) are excluded let us now try 
to understand to what confidence level they can be excluded. The simplest way of exploring 
the parameter space of the magnetized models goes, in short, as follows: 

• the joined two-dimensional marginalized contours for the various cosmological param- 
eters identified already by the analyses of the WMAP 3yr data are ellispses with an 
approximate Gaussian dependence on the confidence level; 

• the confidence intervals for the 2 supplementary parameters of the model (i.e. n^, -Bl) 
can then be determined by using an appropriate gridding approach; 

• the remaining parameters of Eq. (IS. 91) are assumed to be known and follow a Gaussian 
probability densityfunction. 

This approach is rather standard when exploring the impact of new scenarios on the CMB 
observables (see, e.g., [SOl [STl |5l] for the case of non-adiabatic modes supplementing the 
standard ACDM scenario). In the approach we just described, the shape of the likelihood 
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function can be determined by evaluating the least square estimator 



X^(nB,SL) = E 



(obs) 



a 



(obs) 



(6.1) 



where af^^ are the estimated errors from the observations for each multipole and where the 
functional dependence of Ci{n-Q,Bi^ is given by the underlying theory (i.e. the magnetized 
ACDM model) which we try to falsify by comparing its predictions to the actual observations. 
The observed angular power spectra (i.e. Cf^^) are derived by using the (further) estimators 
C7^^ and C7^\ i.e. 



(TT) 



1 



2£ + 



i(T)|2 



c, 



(TE) 



m=-i 



2i- 



\"'lm "'Im I' 



(6.2) 



m=—£ 



whose distribution becomes Gaussian, according to the central limit theorem, when i ^ oo. 
The minimization of Eq. fl6.ll) is equivalent to the minimization of the lognormal likelihood 
function C = —2 InL where L is given by 



L(data| ub, Sl) = L^i^^e 



-x''{«b,Bl)/2 



(6.3) 



Thus, the minimization of Eq. (16.11) is equivalent to the maximization of the likelihood 
of Eq. (16. 3p . In Fig. [11] the parameter space of the model is illustrated, respectively, for 
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Figure 11: The parameter space of the magnetized CMB anisotropics for two illustrative 
choices of the magnetic pivot scale. In both plots /5 = 0. 

two different choices of the magnetic pivot scale /cl, i-e. /cl = Mpc~^ (plot at the left) 
and ki^ = O.lMpc^^ (plot at the right). The shaded spots, in both plots of Fig. [TTl are 
meant to emphasize the value for which the estimator of Eq. (16.11) is minimized; note that 
Xdof — Xmin/-^dof IS the valuc of the (reduced) at the minimum (i.e. when = Xmin)- 
In both plots the data points for the TE and TT angular power spectra have been used 
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in their unbinned form. Overall the total number of data points is A^tt + ^te = 1998 
since A^tt = 999 and N^e = 999. In both plots of Fig. [11] the boundaries of the two 
regions contain 68.3 % and 95.4 % of likelihood as the values for which the has increased, 
respectively, by an amount Ax^ = 2.3 and Ax^ = 6.17. The latter figures stem directly 
from the fact that we are dealing with a two-dimensional parameter space. Figure [TT] offers 
a more quantitative interpretation of the plots reported in Fig. [TUl In Fig. [TD] the value 
of the magnetic pivot scale is /cl = Mpc~^. Therefore the results of Fig. [11] do apply. In 
Fig. [11] the point [ub, -Bl) = (2, 10 nG) is located very far from the outer contour. In a 
frequentist perspective, a model located beyond the outer contours of Figs. [TT] is excluded, 
by the current WMAP data on the TT and TE correlations, to 95% confidence level. 

The exclusion plots reported in Fig. [TT] have been obtained by means of a grid approach 
and by using directly the numerical results. We can now ask ourselves the following question. 
What happens if we use the approximate form of the angular power spectra derived in the 
present paper and use, simultaneously, the same data but in their binned form? The results 
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Figure 12: The parameter space of the magnetized CMB anisotropics as derived by using 
the semi-analytic expressions for the angular power spectra together with the data in their 
binned form. In both plots (3 = 0. 



are illustrated in Fig. [12] It is amusing to notice that, in spite of numerical differences which 
are fully justified the shapes and the regions of the exclusion plots of Fig. [12] are consistent 
with the ones of Fig. [TT] The results of Fig. [12] have been obtained by parametrizing the 
angular power spectra as in Section [5] Note, however, that in Section [5] we only reported 
explicit expressions in the case ric > 1. To obtain the result of Fig. [12] we also need the 
analog formulas but in the casqlj < I- The results of Fig. [T2] show a fair consistency 
of our analytical approach. At the same time it is clear that the value of the reduced 
is larger than in the case of Fig. [TT] in Fig. [T21 the data have been used in their binned 



Recall, as stressed in [W (third paper in the reference), that the cases ub 
separately regularized at the level of the magnetic energy density. 



1 and UB = 5/2 should be 
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form. The largeness of the reduced simply means, within the present approach, that the 
uncertainties entering Eq. (16. ip have been underestimated. With the last proviso, the results 
derived in this paper seem to allow for approximate evaluations of the parameter space of 
the magnetized CMB anisotropies which are compatible with the fully numerical results. 

7 Concluding considerations 

The (limited) question addressed in this paper is how the angular power spectra of the CMB 
anisotropies scale with the parameters of an ambient magnetic field when the remaining 
parameters are close to the ones of the standard ACDM scenario. 

To get a definite answer, the calculation of the magnetized temperature and polarization 
anisotropies had to be reduced to the evaluation of a set of basic integrals whose explicit form 
simplifies in the limit of sufficiently small angular scales, i.e. in the limit of sufficiently large 
multipoles. It has been shown that the temperature and polarization observables obtainable 
by semi-analytic means are sufficiently accurate to infer a set of scaling relations which can 
be used to determine the effect of large-scale magnetic fields on the TT, TE and EE angular 
power spectra. It has been also shown explicitly how the distortions patterns induced by 
large-scale magnetism can be deduced from generalized magnetic form factors accounting 
for the scaling properties of the CMB observables as a function of the parameters of the 
ambient magnetic field. 

The cleanest probe of large-scale magnetism is represented by the polarization autocor- 
relations. Indeed the magnetized EE angular power spectra exhibit a single periodicity, a 
milder dependence upon the dissipative scale an a rather clear shift of the peaks for interme- 
diate and large multipoles. The obtained results have been confronted with the numerical 
calculation. The same evolution equations studied analytically have been integrated numer- 
ically. The analytical derivations are corroborated by the numerical results which are based 
on the same physical description. The obtained results are relevant for dedicated strategies 
of parameter extraction such as the ones mentioned at the end of the previous section. 

Large-scale magnetism is a complicated phenomenon which we choose to scrutinize in 
its infancy, i.e. around the time of photon decoupling. This analysis must necessarily be 
cautious and modest. The rationale for such a caveat does not reside in a particular ethical 
conviction but in the nature of the problem: while in the laboratory we completely dominate 
the initial conditions of our experiments in astrophysics and cosmology this is not the case. 
So, unless we are vigilant it is well possible to mistake an effect due to the peculiar nature 
of initial conditions with a missing piece of dynamics which should have been included in a 
particular regime 0. In this respect it is important to understand in detail all the potentially 

-•^^By this statement I mean that it is important to scrutinize systematically all the relevant plasma effects 
which are applicable in the pre-decoupling regimes. This has been partially done in the present paper as well 
as the original studies of Refs. [HI [3S1 HHl SS] • Even if the present estimates suggest that other (potentially 
relevant) plasma eifects do not play a crucial numerical role, care must betaken in sharpening these estimates. 
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interesting physical effects which could modify the CMB observables. This study is time- 
consuming both at a numerical and at the analytic level. It could even be said, by some, that 
such an approach is slow and pedantic. While we are ready to take this risk, it is rewarding 
that, at the moment, the estimates of the effects of large-scale magnetism on the scalar 
modes of the CMB anisotropies seem to reach (slowly) the same standards employed in the 
absence of ambient magnetic fields. More theoretical effort, in this direction, is certainly 
needed. 



Indeed, a missing piece of the dynamics could be mistaken as a peculiar feature stemming from the initial 
conditions. 
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